The critical role of shear in gravitational collapse 
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We investigate here how the shearing effects present within a collapsing matter cloud influence the 
outcome of gravitational collapse in terms of formation of either a black hole or a naked singularity 
as the final end state. For collapse of practically all physically reasonable matter fields, we prove 
that it would always end up in a black hole if it is either shear-free or has homogeneous density. 
Thus it follows that whenever a naked singularity forms as end product, the collapsing cloud must 
necessarily be shearing with inhomogeneous density. Our consideration brings out the physical 
forces at work which could cause a naked singularity to result as collapse end state, rather than a 
black hole. 
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One of the most intriguing phenomena in black hole 
physics and gravitation theory that has attracted con- 
siderable attention in recent years is the occurrence of 
naked singularities as possible end state of a continual 
gravitational collapse. It turns out that depending on the 
nature of the initial data of the matter distribution and 
the metric functions, cither a black hole or a naked singu- 
larity results as the final outcome of an endless collapse 
(see e.g. yj and references there in for further details). 
In fact, the occurrence of naked singularities could be 
generic because the initial data evolving to a naked sin- 
gularity would form an open set |2] • Such a scenario will 
be of physical interest because a visible extreme strong 
gravity region could possibly provide us with an oppor- 
tunity to observe directly the effects of quantum gravity, 
as generated by the strong curvature regions. 

An obvious question of considerable interest and sig- 
nificance is then, what are the possible physical factors 
which distinguish one outcome from the other. In other 
words, one would like to inquire what are the physical 
agencies operating during the process of a continual col- 
lapse of a massive matter cloud that lead to the forma- 
tion of a naked singularity as the end state rather than 
a black hole, or vice versa. Such an investigation should 
help us understand much better the physics of black holes 
and naked singularity formation in gravitational collapse. 
Towards this end, it was recently shown Q that for the 
spherical dust collapse the shearing effects and inhomo- 
geneity present within the collapsing cloud do play a cru- 
cial role in delaying the formation of the trapped surfaces 
and the apparent horizon. These could in fact cause suf- 
ficient distortion in the geometry of apparent and event 
horizons, thus exposing the singularity to external ob- 
servers. 

Our purpose here is to show that spacetime shear plays 
a crucial role towards determining the black hole or naked 
singularity as end states of gravitational collapse for al- 
most all physically reasonable matter fields. We consider 
here collapse for type I matter fields, which is a rather 
broad class comprising of practically all physically inter- 
esting matter forms such as dust, perfect fiuids, massless 



scalar fields, and such other forms of matter [4|. We show 
that the collapse always ends in a black hole whenever 
there is no shear present in the collapsing cloud. In other 
words, whenever a naked singularity has developed as the 
end state of collapse, it is then always due to the distor- 
tion of the trapped surface geometry caused by the pres- 
ence of shear. As the matter considered here is generic, 
with non- vanishing radial as well as tangential pressures, 
we argue that generically it is the spacetime shear which 
is the physical agency that causes naked singularity as 
the collapse outcome, as opposed to a black hole. 

Consider a spherical collapsing cloud which can be de- 
scribed by the general metric in the comoving coordinates 
{t, r, 9, (f)) as given by, 

ds" = -e^'^'^'^-Ue + e2^(*''^)dr2 + R\t, r)dn^ (1) 

The energy-momentum tensor for any matter field which 
is type I is then given in a diagonal form, 

T/ = -p; T;=p,; T,« = r;=pe (2) 

The quantities p, Pr and pe are the energy density, radial 
and tangential pressures respectively. We take the matter 
field to satisfy the weak energy condition, i.e. the energy 
density measured by any local observer be non-negative, 
and so for any timelike vector V^, we must have, 

T.kV'V'' > (3) 

which amounts to, 

P>0; p + pr>0; p + pe>0 (4) 

It should be noted that we are not confining ourselves 
here to any special type of matter such as e.g. dust or a 
perfect fiuid form, but all forms of matter are included 
where the stress-energy tensor admits one timelike and 
three spacelike eigen vectors. In that sense, our conclu- 
sions are rather generic, applying to a large variety of 
collapse models. 

Now for the metric the Einstein equations take the 
form, in the units (SttG = c = 1) 
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G{t, r) = e-^^{R! f- H{t, r) = e-'"'{Rf (9) 
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In the above, we have introduced the arbitrary function 
F — F{t,r), which has an interpretation of the mass 
function for the cloud, giving the total mass in a shell 
of comoving radius r. We have _F > from the energy 
conditions. 

The shear tensor for the collapsing matter is given by 
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Now let us consider the situation when the matter 
shear vanishes identically. This means. 
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R = qWe" 



(11) 



(12) 



We can use the scaling freedom available in rescaling the 
radial coordinate r, and with a suitable rescaling we can 
always choose q{r) = r. We then have, 

R = re^ (13) 

It then follows that the spacetime geometry (1) becomes. 



^dn^ 



(14) 



Hence we see that there are now five total field equa- 
tions with six unknowns as given by, p, pr, pg, ip, v and 
-F, thus giving us the freedom of choice of one free func- 
tion if we are to complete the solution. Also we require 
the regularity of the initial data, and in particular the 
density distribution, at the initial surface t = ti from 
which the collapse develops. The collapse condition re- 
quires i? < 0, which amounts to ■0 < 0. The singularity 



occurs at R{r,ts{r)) = 0, which implies that in the limit 
as t ^ t g , 



(15) 



It is now possible to integrate the equation (7), using 
equation H13|l . to obtain 



a{t)ij; 



(16) 



where a{t) is an arbitrary function of integration. Since 
the left hand side of the above equation is positive by 
definition, it follows that t/j < implies a(t) < 0. 

Now let us consider the expansion O, for the infalling 
matter congruence. For the metric the expansion is 
given as. 



(17) 



From equation Hll|) and H16|) then it can be easily seen 
that, in case of a shear free collapse. 



e(t,r) = e(i) 
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(18) 



It is obvious that the expansion here is negative, as the 
matter is going through a process of continual collapse. 
As our main interest is to study the curvature singular- 
ity at i? = 0, we require that there are no shell crossings, 
that is, we have R' > in the spacetime. This is equiv- 
alent to the condition that rip' > —1 as we see from the 
expression for R above. We note that the comoving co- 
ordinate system we have employed here is based on the 
congruence of in-falling matter. In this case then ensur- 
ing that there are no shell crossings would ensure that 
the coordinate system is valid and does not breakdown 
till the curvature singularity at i? = 0. 

It is important to note here that in the case of dust, 
which is a special case of general Type I matter fields 
considered here, the congruence of the curves of collaps- 
ing matter are geodesies and hence in that case no shell 
crossings would imply there are no conjugate points in 
the coordinate system. In fact, in this case, we have. 



ip = R'/R' 



(19) 



so i?' ^ implies & ^ —oo, which shows that the shell- 
crosses are equivalent to occurrence of conjugate points in 
the congruence of geodesies. Again the condition R' > 
then ensures that the coordinate system is valid till the 
curvature singularity at i? = 0. In the general case con- 
sidered here, the collapsing matter need not move along 
geodesies, however, imposing the condition as above that 
there are no shell-crossings ensures that the co-ordinate 
system does not break down till the curvature singularity. 

Now for the curvature singulartity at i? = 0, we have 
^ — oo as the singularity acts like a sink for all the 
curves of the collapsing congruence, and the volume ele- 
ments shrink to zero along all the collapsing trajectories 
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0. This implies that at the curvature singularity we 
have a{t) = 0. Further, we note that at all the regular 
spacetime events a{t) is finite and non-zero. This follows 
from equation (11) which implies that ip is finite at all 
regular points because so are R and R. Then from equa- 
tion (16), since vlr, t) has to be regular in the spacetime, 
it follows that a{t) has the above behaviour. This makes 
physical sense also because Q is the physical parameter 
characterizing the volume expansion (or shrinkage) of the 
collapsing cloud. 

Now, let the time taken for the central shell at r = 
to reach the singularity be denoted by tgg , where the sin- 
gularity curve ts(r) corresponds to singularity at i? = 0. 
In the present situation, we need to examine if any non- 
spacelike trajectories could escape away from the singu- 
larity, thus making it visible. In other words, we have to 
find if there are any future directed non-spacelike curves 
which reach to faraway observers, and in the past which 
terminate at the singularity. Suppose now that the sin- 
gularity curve ts{r) is an increasing function with r. Con- 
sider then the spacelikc surface t — tsg- Any event on this 
surface with r > then lies in the spacetime, because the 
singularity R — for this collapsing shell is reached at a 
later time. So this is a regular event at which must be 
finite as shown above. However, this is not possible as 
we know that 6(iso) = — oo. Similar argument applies 
if ts{r) were a decreasing function of r. It follows that 
ts{r) = tgg, which is a constant function, and we have. 



e(i,(r)) = e(t,j--cx) 



(20) 



This implies that the singularity ts(r) is simultaneous. 
The singularity curve being constant is a necessary and 
sufhcient condition for the collapse to be simultaneous. 

A simulataneous collapse, however, necessarily must 
give rise to a covered singularity at i? = 0, and there 
cannot be any outgoing future directed non-spacelike 
geodesies coming out from the same. Because, if there 
were any such outgoing geodesies, given by say t — t(r) 
in the (t,r) plane, which came out from t = ts,r — 0, 
then the time coordinate must increase along these paths. 
This is, however, impossible as there is complete collapse 
at t = tsg, and there is no spacetime beyond that. Hence 
no values t > t^^ are allowed within the spacetime which 
does not extend beyond the singularity. Thus, the col- 
lapse gives rise necessarily to a black hole in the space- 
time (see also 0). Similar argument of course apphes 
either to the points on the singularity curve at r = or 
r > 0. 

It follows that in a shear-free collapse of a general Type 
I matter field, the end state of gravitational collapse is 
always necessarily a black hole. What the above shows 
is that in this case, all the collapsing shells reach the 
singularity simultaneously. As a result, the singularity is 
necessarily covered and is a black hole. 

As we indicated above, we have considered here a suffi- 
ciently general form of matter, and hence our conclusions 
above on the black hole formation as collapse end state 
are rather generic. To get a better insight into how shear 



operates in a dynamically evolving scenario, let us now 
consider in some detail the case of a collapse evolution 
where we assume the matter density to remain homoge- 
neous throughout. We shall construct below an explicit 
class of collapse models to understand how shear works 
and affects the collapse. The choice of a homogeneous 
density profile of course does not mean that we are choos- 
ing any particular form of matter, but it just implies that 
the density is a function of time only, i.e. 



(21) 



Now, let us choose the class of velocity profiles for the 
collapsing shells as determined by the choice. 



v{t,r) = A{R) 



(22) 



Here the function A{R) is any arbitrary, suitably differ- 
entiable function of the physical radius R of the cloud, 
with the initial constraint 

A{R)\t^u = Mr) (23) 
Again, from the Einstein equation Q we get, 

'2(peo-Pro) Pro 



Mr) 



\r{pQ+Pro) PO+Pro 



dr 



(24) 



Let us now assume that the initial pressures have phys- 
ically reasonable behavior at the center r = 0, in that 
the pressure gradients vanish, i.e. p'r^^iO) = p'ogi^) ~ 0, 
and also that the difference between radial and tangential 
pressures vanishes at the center, i.e. Pro(O) ~ PeoiO) = 0, 
which ensures the regularity of the initial data at the cen- 
ter of the cloud. Then, from equation (|24|l . it is evident 
that h'oi^r) has the form. 



t'o(f) = r'^g{r) 



(25) 



where, g{r) is at least a function of r for r = 0, and 
at least a function for r > 0. From equation H25|l we 
can now generalize the form of A(R) as, 

A{R)^R^giiR) (26) 

where gi{R) is a suitably differentiable function and, 

5i(i?)|t=t. =5M (27) 

Now from equations (|21|l and lO we can directly calculate 
the function F and pr as. 



F = ^p{t)R' 



(28) 



and, 



-Pit)-lp(t)j 



Pr ^ -p[t) - -p[t)- (29) 

Also, using equation in equation Q, we have, 

G{t,r) =bir)e^^ (30) 
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Here b{r) is another arbitrary function of the radial co- 
ordinate r. (A comparison with dust collapse models in- 
terprets b{r) as the velocity function for the shells). We 
can write b{r) as, 



6(r) = l + rHQ{r) 



(31) 



Thus we see that for an e ball around the central shell, 
the function G behaves as. 



Now using equation ||SJ) we get 



R 
R 



-e-\l2g,{R) + -p{t) 



(32) 



(33) 



Hence, it is evident that in the vicinity of the singularity, 
that is in the limit i? — > and p oo, and close to the 
central shell. 



R 
R 



= fit) 



(34) 



Here f{t) is another function of time. Thus from equa- 
tion (|34|l and (|29|l we see that in the limit of approach 
to the singularity and near the central shell, the radial 
pressure behaves as. 



Pr = Prit) 



(35) 



Again, in the same limit we can write the tangential pres- 
sure as, 



2pe = RA^ip + Pr) + 2pr « 2pr{t) 



(36) 



It is clear therefore from the above that in the case 
of a homogeneous collapse, a large class of solutions as 
given above, and characterized by the functions A{R) ex- 
ists, for which both the radial and tangential pressures 
homogenize close to the center and in the vicinity of the 
singularity, and hence the collapse becomes necessarily 
shear-free at this limit (note that in spherical symme- 
try homogeneity, through the no energy flux condition, 
implies vanishing shear). The final outcome of such a 
homogeneous collapse is then necessarily black hole. 

We have thus shown that it is black hole whenever the 
collapsing matter is shear-free or homogeneous. That 
is, the end product of collapse could be different from 
black hole (i.e. naked singularity) only if the collapsing 
matter is both shearing as well as inhomogeneous. This 
is an important and a fairly general feature because the 
matter field we have considered is general enough. It 
follows that for a number of classes of general Type I 
matter fields, the homogeneous collapsing configurations 
are subclasses of shear-free collapse. The point we make 
here is a collapse which is homogeneous in density tends 
to a shear-free configuration in the limit of approach to 
the singularity. 

An important insight that emerges from our consider- 
ation here is that there exists a remarkable tying up of 



shear and inhomogeneity for a gravitationally collapsing 
matter cloud to end up in a naked singularity. Intu- 
itively, the physical process that could distort the shape 
of apparent horizon surface so as to expose the singu- 
larity is shear. That is why shear-free collapse always 
ends in black hole. Further we have also shown that 
as the singularity is approached, shear could only be 
produced by inhomogeneity in density. Note that this 
is also the case for pressureless dust (and homogeneous 
pressure) where shear is solely generated by density in- 
homogeneity. What seems to happen is that as the sin- 
gularity is approached pressures turn homogeneous and 
hence ineffective. Thus inhomogeneous density is neces- 
sarily required for shear to be non-zero as the singularity 
is approached. Interestingly, the same tying up of inho- 
mogeneity and shear (anisotropy) occurs for non-singular 
cosmological models as was first argued by one of us P] . 
We hence see that the collapse for a general matter field 
as we have considered here generically tends to the out- 
come which is observed in the dust case Q. This raises 
an interesting possibility that generically all collapse con- 
figurations could tend to a dustlike model in the vicinity 
of the singularity. 

Another interesting feature that seems to follow quite 
generally and generically (where weak energy condition 
is satisfied with densities and pressures being positive), 
when one considers the outgoing null geodesies equation, 
is that as expected it is only the central singularity at 
r — Q which could be naked and none else for r > 0. 
To see this, and in order to examine the nature of the 
central singularity ai R — Q,r — consider the equation 
for outgoing radial null geodesies. 



dt 
dr 



(37) 




One could now write the above equation in terms of the 
variables [u — r^jR), where one could choose a = |, 
and using equation we get, 



(38) 



where we have written R = rv. Now if the null geodesies 
do terminate at the singularity in the past with a definite 
tangent, then at the singularity we have ^ > 0, in the 
(m, R) plane with a finite value. It follows that all points 
r > on the singularity curve are then covered because 
F/ R — > oo (when weak energy condition is satisfied and 
pressures are positive then F(r) tends to a finite posi- 
tive value for any r > on the singularity curve) with 
^ ^ — oo. This implies that non-central singularities 
will always be covered for a general enough type /matter 
distribution. The central singularity could however be 
visible. 

More generally, it can be stated that any spacetime 
or evolving distribution has two dynamically interesting 
features, namely the homogeneity and isotropy. Though 
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they are quite different properties yet they seem to get 
quite intimately correlated at the end stage of the col- 
lapse in the vicinity of the singularity. That is, as the 
singularity is approached one implies the other. The role 
of shear or anisotropy is to make the motion incoherent. 
This means the motion is not in the same sense in all 
directions, and the collapse or expansion may proceed 
at different rates in different directions. Intuitively, the 
shearing force will act against the concentration of mat- 
ter quickly enough in small enough region, or creates a 
distortion, which in turn delays the formation of apparent 
horizon or the trapped surfaces. This is how the shear, 
which in the ultimate stage of collapse entirely hinges 
on density inhomogeneity, plays the key role. Thus both 



shear and inhomogeneous density are required if the col- 
lapse has to end in a naked singularity. 

What we have shown here is that for a fairly generic 
matter distribution, the absence of shear would always 
lead to a black hole. Hence in a way we have pointed out 
and identified here the physical mechanism that causes 
naked singularities in collapse. In this sense, our results 
here explain in a natural manner why naked singularities 
do develop as end state of collapse, when there are 
non-vanishing shearing forces and inhomogeneities 
present within the collapsing cloud. 

It is a pleasure to thank Sanjay Jhingan and Louis 
Witten for discussions. 



[1] G. Horowitz, |hep-th/0312T23| A. Krolak, Prog. Theor. 

Phys. Suppl. 136, 45 (1999); R. Penrose, in Black holes 

and relativistic stars, ed. R. M. Wald (University of 

Chicago Press, 1998); R. M. Wald, gr-qc/9710068 P. S. 

Joshi, Pramana 55, 529 (2000); M. Celerier and P. Szek- 

eres, gr-qc/0203094 
[2] T. Hertog, G. T. Horowitz and K. Maeda, gr-qc/0307102" 

P. S. Joshi and I. H. Dwivedi, Class. Quantum Grav. 16, 

41 (1999). 

[3] P. S. Joshi, N. Dadhich and R. Maartens, Phys. Rev D65, 

p. 101501R (2000). 
[4] S. W. Hawking and G. F. R. Ellis, The large scale structure 

of spacetime, Cambridge Univ. Press, Cambridge (1973). 
[5] D. Kramer, H. Stephani, E. Herlt, and M. MacCallum, 



Exact solutions of Einstein's field equations, CUP, Cam- 
bridge, (1980); for a further discussion on shear in the 
context of gravitational collapse, see also E. N. Glass and 
B. Mashhoon, Ap. J. 205, p. 570 (1976); R. A. Sussman, J. 
Math. Phys. 29, p 945 ((1988); E. N. Glass, J. Math. Phys. 
20, 1508 (1979); E. Brinis, S. Jhingan and G. Magli, Class. 
Quant. Grav. 17, 4481 (2000); S. M. C. V. Goncalves, 
gr-qc/0312123 

[6] N. Dadhich, J. Astrophys. Astr. 18, 343 (1997); N. Dad- 
hich and A. K. Raychaudhuri, Mod. Phys. Lett. A14, 2135 
(1999). 

[7] R. Goswami, P. S. Joshi, gr-qc/0310122l To appear m 
Phys.Rev.D 



